Atom-ion quantum gate 
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We study ultracold collisions of ions with neutral atoms in traps. Recently, ultracold atom-ion 
systems are becoming available in experimental setups, where their quantum states can be coherently 
controlled. This allows for an implementation of quantum information processing combining the 
advantages of charged and neutral particles. The state-dependent dynamics that is a necessary 
ingredient for quantum computation schemes is provided in this case by the short-range interaction 
forces depending on hyperfine states of both particles. 

In this work we develop a theoretical description of spin-state-dependent trapped atom-ion col- 
lisions in the framework of a Multichannel Quantum Defect Theory (MQDT) and formulate an 
effective single channel model that reduces the complexity of the problem. 

Based on this description we simulate a two-qubit phase gate between a '^^^Ba^ ion and a *^Rb 
atom using a realistic combination of the singlet and triplet scattering lengths. We optimize and 
accelerate the gate process with the help of optimal control techniques. Our result is a gate fidelity 
1 - 10"^ within 350iis . 

PACS numbers: 34.50.cx,37.90.-|-j,03.67.Bg,03.67.Lx 



I. INTRODUCTION 

Ongoing developments in quantum information pro- 
cessing stimulate an intense search for physical systems 
suitable for its implementation. Beside solid-state and 
photonic systems, cold ions and neutral atoms represent 
major candidates in this direction. 

Neutral atoms can be accurately manipulated in dipole 
traps P, 01, optical lattices 3 or on atom chips . 
Advanced evaporative and laser-cooling techniques allow 
their preparation in the vibrational ground state of dif- 
ferent trapping potentials. Single ions can be confined 
in Paul or Penning traps @ and sideband laser cooling 
allows to cool them down to the trap ground state. 

Ultracold systems combining ions and neutral atom are 
currently being explored 0, H, Hi- Besides several new 
quantum mechanical aspects of this system, the studies 
are motivated by potential applications. For example 
techniques of sympathetic cooling of trapped atoms by 
laser-cooled trapped ions can be developed [l^, [ll| . In 
this paper we propose a scheme for a quantum gate that 
combines the advantages of atoms and ions for quantum 
computation. 

The trapping potentials of atoms and ions, although 
both built up with oscillating electromagnetic fields, do 
not interfere with each other, since the oscillation fre- 
quencies of the respective fields typically differ by orders 
of magnitude. The strength of the effective ion poten- 
tial can be much stronger than for neutral atoms. Tight 
confinement enables fast transport and together with the 
good addressabihty of single trapped ions with lasers, this 
is among the advantages of ions for implementing quan- 
tum computation. 

Realization of the Mott insulator phase allows to pre- 



pare an array of atoms with well controlled number of 
particles in a single site of an optical lattice. In compar- 
ison to ions 

The possibility to prepare an array of atoms in the 
Mott insulator phase in an optical lattice combined with 
the long decoherence times of neutral atoms is a reason 
to use atoms for the storage of quantum information. 
Furthermore, the two-particle interaction of an atom and 
an ion is typically much stronger than for two neutral 
atoms, which allows fast gate operations. 

While qubits can be stored in internal electronic de- 
grees of freedom of both kinds of particles, the state de- 
pendent dynamics suitable for two-qubit gates requires 
engineering of the two particle interaction. To this end 
one can use external electromagnetic fields, e.g. magnetic 
Feshbach resonances that allow for a precise tuning of the 
two-body effective scattering properties. The long-range 
atom-ion interaction also supports a trap-induced type 
of resonance [l^l , because of their generally state depen- 
dent nature, they will constitute a basic element of our 
quantum computation scheme. 

In this work we solely make use of these trap-induced 
shape resonances that occur at relatively large distances. 
In this way we avoid some possible unwanted processes 
that may result from molecular dynamics at short dis- 
tances. We nevertheless plan to include magnetic Fesh- 
bach resonances to our theory, which can be applied to 
perform two-qubit operations in a controlled collisions 

A possible setup for quantum computation is schemati- 
cally depicted in Figure[T] Atoms are stored in an optical 
lattice in a Mott insulator phase such that each lattice 
site is occupied by exactly one atom. One movable ion is 
used to create long-distance entanglement between pairs 
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of atoms and to perform quantum gates. The basic in- 
gredient of this idea is the controUed and qubit-sensitive 
interaction between atoms and ions. In this paper wc 
focus on the dynamics of a single atom interacting with 
a single ion, nevertheless; our approach can be easily ex- 
tended to the situation of many atoms, or more than one 
ion, at a later stage. 




atoms in optical lattice 



FIG. 1: (Color online) Concept for quantum computation 
with atoms and ions: Atoms are prepared in an optical 
lattice in a Mott insulator phase. A movable ion en- 
tangles the atoms and can also be used for sympathetic 
cooling. 

In this work we develop a theoretical model for spin- 
dependent atom- ion collisions. For the case of an 
alkaline-earth ion and an alkali atom wc formulate a 
model based on the Multichannel Quantum Defect The- 
ory (MQDT) 0] taking into account the presence of trap- 
ping potentials Within our model the atom-ion in- 
teraction is described by the long-range 1 /r^ polarization 
potential combined with a set of quantum-defect param- 
eters representing the effect of the short-range potential. 
The essential parameters for our approach are the singlet 
and triplet scattering lengths, which arc not yet known 
with sufficient accuracy, but can probably be measured 
in upcoming experiments. In our paper we discuss sev- 
eral different regimes for the values of the singlet and 
the triplet scattering length. For some specific range of 
scattering lengths we are able to reduce the complexity 
of the problem by employing an effective single-channel 
model for the atom-ion dynamics. 

Using the effective single-channel description we are 
able to simulate a two-qubit phase gate for an arbitrary 
combination of atom and ion species. Applicability of 
the model, however, requires values of singlet and triplet 
scattering lengths that are nearly equal. In this case, 
within the single-channel description and for a specific 
system of ^^'^Ba"'' ion and ^^Rb atom, we develop a phase 
gate process yielding a fidelity of 1 — 10"'^ within the 
gate time of 346 ps. Being equivalent to the CNOT 



gate, the phase gate is universal for quantum compu- 
tation [l^. Therefore we demonstrate the feasibility of 
quantum computation on the system under considera- 
tion. In a general situation, when the scattering lengths 
are not similar, the phase gate can be even faster; how- 
ever this requires going beyond the single-channel effec- 
tive description and is outside the scope of the present 
paper. 

There are two main mechanisms that could lead to a 
failure of the quantum gate. One is the radiative charge 
transfer, which in our scheme leads to a loss of both par- 
ticles in case of heteronuclear species. In contrast, for a 
homonuclear collision 0] the charge transfer results in a 
physically equivalent situation and therefore cannot be 
considered as a loss mechanism. Heteronuclear alkaline- 
earth ion - alkali atom systems have the advantage of 
a relatively simple electronic level structure, and for the 
systems studied so far; Na-Ca"*" [§, [TH and Rb-Ba+ [l7| . 
the charge exchange rate remains much smaller than the 
elastic coUisional rate, even in the presence of resonances. 
The second loss process results from spin changing col- 
lisions. In our scheme the qubits are encoded in hy- 
perfine spin states, and collisions leading to final states 
outside of the computational basis have to be avoided. 
In the regime of applicability of our single-channel ef- 
fective model, the coupling between different channels is 
by definition very weak and those kinds of losses can be 
safely neglected. Even in a general situation, anyway, a 
multichannel treatment including all possible spin-state 
channels offers the possibility to gain control over spin- 
changing processes by appropriate engineering of the gate 
dynamics. 

The paper is organized as follows. In Sec |TT] we de- 
scribe the basic setup and model used throughout the 
paper. Further we briefly discuss the atom-ion polariza- 
tion interaction and we introduce the concepts of cor- 
relation diagrams and trap induced resonances. The 
MQDT for trapped particles as well as its reduction to 
a single-channel model is developed in Sec IIIII The pre- 
sented theory allows for computation of eigenstates and 
eigencnergies either in single-channel or multi-channel 
situations. The dynamics of atom-ion collisions is dis- 
cussed in Sec. IIVI where correlation diagrams and the 
Landau Zener theory help to understand the features of 
the system. Sec|V] presents the concepts and results of 
our two-qubit phase gate simulations. A summary of our 
results, together with further perspectives and ideas, is 
given in Sec IVII 



II. BASIC SETUP AND MODEL 

We consider a system consisting of a single atom and 
a single ion, stored in their respective trapping poten- 
tials. Such potentials can be created with rapidly oscil- 
lating (rf) electric fields for ions and with optical traps 
based on the ac Stark effect for atoms. These traps can 
be well approximated as effective time-independent har- 
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monic traps, as long as the particles are close to the 
ground state of the potential. For this setup we intro- 
duce an effective Hamiltonian 
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where ^^(a) is the mass of the ion (atom), tona) a-nd di(^a) 
denote frequency and location of the harmonic trapping 
potential of the ion (atom) and W{r) is the interaction 
potential. A microscopic derivation of the Hamiltonian 
Eq. ([T]) can be found in [l^. Here, for simplicity we 
have assumed spherically symmetric trapping potentials 
and the same trapping frequencies for atom and ion: 
uji = tUa = to. We stress, however, that our approach can 
be easily generalized to anisotropic trapping potentials 
and different trapping frequencies |12| . A general treat- 
ment would imply coupled ccnter-of-mass (CM) and rela- 
tive degrees of freedom, thus a six dimensional equation, 
but there are no fundamental difficulties. In experiments 
both traps can be designed to be spherically symmetric, 
while the assumption of the same trapping frequencies al- 
lows us to decouple the relative and CM motions, thereby 
reducing the dimensionality of the problem from six to 
three. This choice simplifies our numerical calculations, 
on the other hand it allows to capture the most important 
features of the system. 

We transform the Hamiltonian Eq. ([T]), introduc- 
ing CM and relative coordinates, Rcm = {rriiri + 
iTiara)/{mi + nia) and r = — r^, respectively. With- 
out loosing generality we can choose the coordinate frame 
such that the vector of trap separation d = d^ — = de^ 
points in the z direction. In this way we obtain the rela- 
tive Hamiltonian 



(0) 



-fiUJ^cf — fiuj'^dz, (2) 



where fi = mima/(jni + nia) denotes the reduced mass 
of the atom-ion system and 



hI"J = -^^r + ^f^^'^r' + Wir) 
is the Hamiltonian for the special case d = 0. 
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Atom-Ion System 


R*{ao) 


lo{ao) 


E* /hikUz) 


1353^+ ^ 87-^^ 


5544 


826 


1.111 


«Ca+ + «^Rb 


3989 


1178 


4.142 


«Ca+ + 23Na 


2081 


1572 


28.545 



TABLE I: Characteristic length and energy scale for 
example systems. Oscillator lengths are calculated with 
uji = uja — in X lOOkHz. 



interaction constant C4 can be expressed in terms of the 
electric dipole polarizability a of the atom in the elec- 
tronic ground state (S-state): C4 = ae'^/2. The electron 
charge is denoted as e. At short distances the interaction 
is dominated by the exchange forces, and higher order 
dispersion terms (Cg/r^, Cg/r®) also become relevant. 
In our approach we model the short-range part of the 
potential using the quantum-defect method, that is we 
do not require the knowledge of the exact form of the 
short-range interaction. The interaction potential, in ad- 
dition to the model potential, is depicted schematically 
in Fig. H 



Ion 



Wir) 




Polarized atom 
— R* 



Vir) = -C4/T* 



FIG. 2: (Color online) The long-range part atom-ion in- 
teraction potential equals —C^/r'^. At distances smaller 
than the potential minimum r^nin, repulsive terms start to 
dominate. Quantum defect theory replaces the actual po- 
tential W{r) (solid line) with a reference potential V{r) 
(dashed line) and includes the short-range effects using 
a quantum defect parameter related to the short range 
phase of the relative wavefunction. The characteristic 
range R* of the interaction is typically much larger than 



Atom-ion interaction 



At large distances the atom-ion interaction potential 
has the asymptotic behavior W{r) ~ — C4/r^ (r —> 00). 
This results from the fact that the ion charge polarizes 
the electron cloud of the atom, and the induced dipole 
and the ion attract each other. Therefore, the atom-ion 
interaction falls into an intermediate category, between 
the long-range Coulomb forces W{r) ^ 1/r and the van 
der Waals forces W{r) ^ for neutral atoms. The 



By equating the interaction potential C^/R** to the ki- 
netic energy we can define some characteristic 
range R* = \j2\iC^IW- and corresponding characteristic 
energy E* = /i^/2/xi?*^ of the atom-ion interaction. Ta- 
ble [I] gives the characteristic range and energy for some 
example atom-ion systems. For comparison it also in- 
cludes the harmonic oscillator length Iq = y^h/ fiuj for 

= Wq = 27r X lOOkHz. 
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B. Single-channel quantum defect treatment 

The short range interaction potential between atom 
and ion is typicaUy quite complicated and in most cases 
it is not known theoretically with an accuracy sufficient 
to determine the scattering properties in the limit of ul- 
tracold energies. In order to avoid complications while 
using the explicit form of the short-range potentials, we 
resort to the quantum-defect method, which allows to in- 
clude the effects of the short-range forces in an effective 
way. This consists in substituting the actual potential by 
the reference potential V{r) = —d/r'^ at all distances 
(see Fig. [5]), and assigning an appropriate short-range 
phase to the wavefunction to model the effects of the 
short-range potential. 

We illustrate this approach by solving the relative 
Schrodinger equation H^'^^^ir) = £;*(r) at d = 0. To 
this end we apply the partial wave decomposition 

^{r)=Y,yirnir)Mr)/r, (4) 

Im 



to obtain the radial Schrodinger equation 



2/.t 9r2 2/i 



^i(r) =0 



(5) 

for the radial wavefunctions Here, are the 

spherical harmonic functions describing the angular part 
of the 3D wavefunction, where I and m are the quan- 
tum numbers of the relative angular momentum and its 
projection on the symmetry axis z, respectively. In the 
limit of r — > wc can neglect trapping potential, energy 
and centrifugal barrier in comparison to V{r) ~ —C^jr^ ^ 
which yields 



2/i 9r- 



with the solution 



C4 
^4 



V,'/(r) = 0, 







(6) 



(7) 



where is a parameter that can be interpreted as the 
short range phase. In our method we treat Eq. ([7]) as a 
boundary condition that we impose on the radial wave 
functions at short distances, while solving the relative 
Schrodinger equation in the general case d 7^ 0. 

In the absence of a trapping potential, and for I = 
and £^ = 0, the solution Eq. ([7]) becomes valid at all 
distances. By comparing the long range behavior of 
Eq. 0: lim^^oo '00(f)/'' ~ {\ ^ R* j r cot ip) , with the 
well known asymptotic form of the s-wave radial wave 
function: linir^oo V'o (?")/'' ^ 1 — a/r, we can relate the 
short range phase to the scattering length 

a — — i?* cot </?, 



(8) 



which is a measurable physical quantity. We note that 
R* determines the typical length-scale for the scattering 



length. In this chapter we have focused only on single- 
channel collisions, not considering internal states of the 
particles. The present approach is generalized in Sec. IIIII 
to the realistic multichannel situation. 



Correlation diagram and trap induced 
resonances 




FIG. 3: (Color online) An example correlation diagram 
calculated for R* = 3.68^o and the short-range phase 
if = 0.747r showing the energy spectrum versus the trap 
separation. The partial wave number / is given for the 
lowest states at = 0. As explained in the text, molec- 
ular state energies have an approximate parabolic d- 
dependence indicated by the dashed line. 

In order to obtain an intuitive understanding of atom- 
ion collisions, we describe them in terms of correlation 
diagrams, showing the energy spectrum as a function 
of the trap separation d (see Fig. [J). Such correla- 
tion diagrams in our case connect the asymptotic vi- 
brational states for large trap separation to the molec- 
ular and vibrational states at zero trap separation. At 
large distances we find harmonic oscillator-like equidis- 
tant eigenenergies that are independent of d. Molec- 
ular bound states, that corresponds to the eigcnstates 
with energies well below the zero-point vibration energy 
Eq = 3/2hLu, experience a quadratic shift with the dis- 
tance d. This can be easily understood by noting that 
the bound states 'I'moi('') are well localized around r = 0, 



and (fmoili^if l^'moi) ~ £'moi + ^fJ-uj'^d'^, where E,-noi is 
the molecular binding energy at d = 0. Beside the given 
arguments the quadratic shift in the molecular energy 
becomes immediately clear in Fig. 21 The molecular po- 
tential 'hangs' from the relative trapping potential in the 
low-distance region, thus increasing (decreasing) d shifts 
up (down) the molecular energy as d^. 

At some particular distances, the energies of the molec- 
ular states become equal to the energies of the vibrational 
levels (see Fig. S]), and the spectrum exhibits avoided 
crossings, known as the trap-induced shape resonances 
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|18l |. By slowly changing the trap separation d we can 
pass through the resonance adiabatically, converting the 
trap vibrational states into molecular states, thus pro- 
ducing molecular ions. Since this process is reversible, 
we can coherently control the dynamics of our system by 
appropriately adjusting the trap distance. 




FIG. 4: (Color online) Trap-induced shape resonance: at 
a certain trap separation d = d^cs , the energy of a molecu- 
lar bound state becomes degenerate to a trap vibrational 
energy. The adiabatic eigenenergies exhibit an avoided 
crossing at this position. The arrows indicate that the 
molecular energy is shifted, if the relative trap position 
is changed. 



III. QUANTUM-DEFECT THEORY FOR 
TRAPPED PARTICLES 



A. Multichannel formalism 

In general, the interaction properties depend on the 
internal state of two colliding particles. For these in- 
ternal states we choose a convenient basis, in which the 
two-particle Hamiltonian is diagonal at large particle dis- 
tance, where the interaction potential is negligible. We 
then refer to the two-particle basis states as scattering 
channels. The wavcf unction is decomposed into the cho- 
sen basis, which allows us to write down Schrodingcr's 
equation in matrix form. In the following we introduce 
an MQDT formalism, following closely the formulation 
by F. Mies [l^, and adopting it to a situation includ- 
ing an external trapping potential. Assuming the same 
trapping frequencies for atom and ion, the CM and rel- 
ative degrees of freedom are decoupled. In this case we 
can describe the relative motion with the close-coupled 
Schrodingcr equation 

- ^A*(r) + [W(r) + U(r) - El] *(r) = 0. (9) 

Here I denotes the identity matrix, W(r) is the interac- 
tion matrix, which is asymptotically diagonal 



(10) 



with indicating the channels. The trapping poten- 

tial U(r) is diagonal at all distances 



;7y(r) = i^w2(r-d)25y. 



(11) 



The matrix '4'(r) contains linearly independent solu- 
tions, where N is the number of channels. The threshold 
energies for the molecular dissociation in channel i are 
denoted by E°°. 



1. Special case: d — Q 

For d = the external potential is spherically sym- 
metric and the dynamics for different relative angular 
momenta / is decoupled. We can decompose F(r) into a 
partial wave expansion 



*(r)=^QFi(r)rro(f)A, 



(12) 



where Here, for simplicity we consider only the m 
subspace. The radial wave functions F; (r) fulfill 



with 



2fi dr- 



2^tr2 ^ 2 



2 2 



(13) 



(14) 



In our calculations we model the short-range potential 
by choosing appropriate short-range phases ipi for each 
of the channels. This is equivalent to setting Wij{r) = 
{E°° — C4/r*) 6ij at all distances. In this case the refer- 
ence potentials that are necessary to define the MQDT 
functions can be simply taken as diagonal elements of 



W(r): VS,r) = 
potentials, Vi{l,r 



Given the reference 



one can associate to them a pair of 
linearly independent solutions fi{l,r) and gi{l,r) of the 
single-channel Schrodingcr equation that have WKB-like 
normalization at small distances 



Ul,r) = fc,(?,r)-i/2sinA(Z,r), 
g,{l,r) ^ h{l,r)-^/^ cos (3,{l,r), 



r ~ rmin, (15) 



where ki{l, r) = ^2[i(E — Vi(l,r)/h is the local wavevec- 
tor and j3i{l,r) = J^dx ki{l,x) is the WKB phase. Here, 
^min denotes a typical distance where the minima of the 
realistic potential occur, and the semiclassical approxi- 
mation is applicable. In our modeling rmin 0, and 
Eq. (fT5|) describe the asymptotic behavior r 0. 

The solution to Eq. (fT3|) can be expressed in terms 
of a pair of functions fi(r) = {Sijfi{l,r)} and gi{r) = 



Fiir) = {i{r)+giir)YiiE) 



(16) 
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Here, Yi{E) is the quantum-defect matrix that repre- 
sents the effects of the short-range potential, in partic- 
ular couplings between channels, and will be discussed 
later. The Matrix A has constant coefficients and is de- 
termined by the boundary conditions at r — > oo. We note 
that in MQDT the functions f;(r) and g;(r) describe in 
general only the asymptotic (r —^ oo) behavior of Fi(r). 
Due to our choice of W(r) and Vi{l,r), however, in our 
case these functions will be valid at all distances. 

In analogy to MQDT in free space, we introduce an- 
other type of solutions, normalized at r oo. At large 
distances the harmonic potential dominates and the so- 
lution vanishing at r — > oo reads 



(17) 



where D^{z) 



is the parabolic cylinder function, E = 
^), and ^ is the harmonic oscillator length 
The two types of solutions, Eq. (fT5|) and 
Eq. ini), can be related by the MQDT functions Vi{l,E) 
aRdMi{l,E): 



(18) 



0,(/, r) = MS, E)ycosv,{l, E)Ml, r) 

- siniy,{l, E)gi{l,r) 



The function Vi{l,E) mixes the two solutions Eq. (|T5l) 
leading to the exponentially decaying function (j)i{l,r), 
whereas Afi{l,E) provides the overall normalization. In 
fact the normalization can be calculated directly from 



MS,E) 



dvi{l,E) 
2n dE 



-1/2 



(19) 



Now, the wave function F;(r) can be equivalently ex- 
pressed in terms of solutions ^i{r) = {5ij4>\{ry\ normal- 
ized at infinity: , 



¥i{r) = ^i{r)A. 



(20) 



By comparing Eq. ^6|) with Eq. ((20|) one arrives at the 
following equation: 

[Yi (E) + tan i^i (E)] Ni (E) cos {E)A = 0, (21) 

where i^i{E) = {S.^iyHE)} and Ni{E) = {S,,MUE)}. 
This has a nontrivial solution (A 7^ 0), if 

\Yi{E)+t&ni^i{E)\^0, (22) 

which is a standard condition determining bound states 
in the MQDT approach. From Eq. ([2^ one can evaluate 
eigenenergies in the multichannel case, while the cigcn- 
states arc given by Eq. (I20p . with A determined from 
Eq. (PT|) . This procedure yields a set of eigenfunctions 
Fin{r) = ^i{r)hin and corresponding eigenenergies Ein, 
where bi„ is a constant vector, and the label n enumer- 
ates the solutions: 



(23) 



Similarly to ordinary scalar wave functions, the mul- 
tichannel eigenstates corresponding to different non- 
degenerate eigenenergies are orthonormal: 

poo 

/ drFUr)^Fi„,{r)=6„,n- (24) 
Jo 



2. Generalization to d ^ 

At nonzero trap separation, the Hamiltonian is no 
longer rotationally invariant, and the procedure pre- 
sented in the previous section based on decoupling of 
states with different values of / does not apply. Never- 
theless, we can utilize the previous solutions at d = to 
diagonalize the full problem Eq. ([9|) at d ^ 0. To this 
end the total wave function is decomposed in terms of 
arbitrary expansion coefficients q„: 



F(r) =^Q„F^„(r)l^o(f)A■ 



(25) 



In 



Substituting this into Eq. ([9]) and setting d = de^ we 
arrive at the set of coupled equations 

{Ei„ + ^fiuj^d^) cin H- ^J^u''dY^ Dfci'n' = Ecln, (26) 

I'n' 

which in principle can be solved numerically with stan- 
dard methods for matrix diagonalization. Here 



Dif = {Yio\ cos 9\Yi,o) drFz„(r)VFz„,(r), (27) 
Jo 



and 



(Viol cos 9\Yi,o) = ldnYilir)coseYi,o{r) 
l + l 



^{2l + l){2l + 3) 
I 

^ V(2^-i)(2zTTy 



(28) 



3. Parametrization of Y; {E) and the frame transformation 

In the regime of ultracold collisions, the variation of 
the total energy E and the height of the angular mo- 
mentum barrier (for the lowest partial waves which are 
important in the ultracold regime Q) are much smaller 
than the depth of the potential at r ~ r^^^ where the ma- 
trix Yi{E) is defined. Therefore it is justified to neglect 
the dependence of Yi{E) on both energy and angular 
momentum, and to set Yi{E) = Y. In this way, deter- 
mining the matrix Y at a single value of energy, we may 
describe the atom-ion collisions in the whole regime of 
ultracold temperatures. 

In the paper we consider collisions of an alkali atom 
with an alkaline-earth ion in their electronic ground 
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states. Hence the asymptotic channel states that are 
used in the Schrodinger equation Eq. p3p . can be char- 
acterized by the hyperfine quantum numbers fi ,m and 
/2,m/2 for ion and atom respectively, and by the angular- 
momentum quantum numbers I and m; of the relative 
motion of the atom and ion CM. In the rest of this sec- 
tion we label those channels by a = {/i/2™/iW/2^"^i}- 

At short distances, the potential matrix becomes di- 
agonal in the molecular basis characterized by the total 
electron and nuclear spins and their projections, because 
the short-range forces depend on the electronic configu- 
ration of the entire atom-ion molecular complex. In fact, 
the molecular potentials that correlate with atom and ion 
electronic ground states at large distances depend only on 
the total electron spin S [§| . For our choice of species the 
electronic configurations are identical as in collisions of 
two hydrogen atoms. Thus, S can take the values (sin- 
glet configuration) and 1 (triplet configuration). Hence 
the quantum defect matrix Y, which contains the full 
interaction information, can be parameterized with only 
two constants, the singlet and triplet at scattering 
lengths. These constants depend only on the species. 

In our approach we apply a frame transformation to 
find Y in the basis of hyperfine states [l^, [2l| . As shown 
in Ref . Q , this approximation is very accurate for atom- 
ion collisions due to a clear separation of length scales as- 
sociated with the short-range and long-range forces. On 
the one hand, exchange interaction becomes significant 
only at distances of the order of few tens of ag (atomic 
units), when the electronic wavefunctions of atom and 
ion begin to overlap. On the other hand, the polariza- 
tion forces are very long-ranged and they are modified 
by the presences of the centrifugal barrier only at large 
distances of the order of R*. 



B. Reduction to an effective single-channel model 

in the case of as ~ at 

The off-diagonal matrix elements in the quantum de- 
fect matrix Y are proportional to the "coupling" scat- 
tering length, that is defined as l/oc = l/os — I /at @. 
Therefore in the case of « at the coupling between 
channels is weak and the many channel description can 
be effectively reduced to a single-channel model. To this 
end we solve the multi-channel problem at d = 0, and 
we find the corresponding eigenenergies Eni and eigen- 
states F„((r) from Eq. ^ and Eq. (QH). If the mixing 
between channels is weak, in each of the eigenfunction 
F„i(r), there is only one channel that dominates, i.e. 
the vector b;„ has only one element which is close to 
unity. We divide the total multichannel spectrum into N 
distinct subsets according to the channel that gives the 
leading contribution to F„;(r), and for each of the sub- 
set we determine the effective short-range phase (^°*^ (or 
the scattering length a^^ = — i?* cot (p°*). This is done 
by matching the multichannel spectrum Eni in each of 
the subsets with the single-channel generated by Eq. ^ 



with the quantum-defect parameter ipf . In the limit of 
zero coupling between channels, the effective scattering 
length a®*^ is equal to a,, = «<. In the presence of weak 
coupling, is in general different from as and at , since 
the asymptotic channels typically correlate both to sin- 
glet and triplet molecular states at small distances. This 
procedure yields a set of N short range phases {(^sj^^}, 
which are used at a later stage in the single-channel calcu- 
lations. We note that the effective phases {'/Sj'*^} depend 
on the relative angular momentum, and in principle they 
weakly depend on the energy. We have verified, how- 
ever, that within the considered range of energies limited 
to the bound states close to the dissociation threshold, 
and to a few tens of the lowest vibrational states, the 
variations of ipf^i^E) with the energy are negligible. 

For collisions when only a single open channel exists, 
the remaining closed channels are typically only weakly 
coupled to the open channel (apart from the case of 
resonances), and the resulting multichannel wave func- 
tion is dominated by the open-channel component. The 
situation changes, however, when there are more open 
channels, and the channel mixing can be significant. We 
have investigated this issue numerically, picking the spe- 
cific ion and atom pair ^"^^Ba^-^^Rb and the trapping 
frequency = Wa = 27r x 30kHz. We have consid- 
ered collisions within the mp = Z subspace, assuming 
that initially the particles are prepared in the channel 
ai = {fi — l,?Ti/i = 1,/a = 2, mj^ = 2}. This choice 
is relevant for our modeling of the quantum gate, as we 
show later. For the collision energies above the dissocia- 
tion threshold of the channel ai, a second open channel 
exists, with a2 = {fi = = 2, fa = l,mf^ = 1}. 

In this case we find that admixture of the two remaining 
closed channels is negligible, whereas the contribution of 
both channels ai and a2 in the multichannel eigcnstates 
is typically large, and the contributions from ai and Q!2 
cannot be separated. The only exception is the case of 
similar as and at, where the inter-channel coupling l/oc 
is small, and the multichannel eigenstates are dominated 
by the single-channel contributions. 

We estimate the validity of the single channel approxi- 
mation by calculating the overlap of the exact multichan- 
nel and the single channel wave functions, in the range 
of trap separations d, interesting for our dynamics. The 
minimum over d yields some overall fidelity related to the 
reduction to the single channel model. For similar singlet 
and triplet phases, we have (fif^ ^ (j)s ^ 'Pt, and the mul- 
tichannel wave functions differ from their single-channel 
counterparts by the presence of negligible contributions 
in the channels others than the dominating one. The 
relative contributions of individual channels are given by 
the vectors bi„ that are obtained in the calculation of 
eigenstates at d = 0. When d > one has to take into 
account that the multichannel wave functions are linear 
combinations of the solutions at d = (see Eq. ([25])), 
and the overlap between single- and multi-channel eigen- 
states is a linear combination of the overlaps calculated 
at d = with the expansion coefficients given by |q„P. 
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Resonances occurring at d = can lead to significant 
channel mixing for few states, although the average fi- 
delity is very high. The more these highly mixed states 
contribute to the wavefunctions, the larger is the error of 
the model. For the error estimation used in this work, 
we only use the trap ground state at maximal d and a 
molecular state at minimal d of interest. We find that the 
two fidelities lie in the same range, thus we take the min- 
imum of them and assume the result as a lower bound 
for the fidelity at intermediate distances. We have ad- 
ditionally verified that reductions to the single channel 
model works best for positive values of singlet and triplet 
scattering lengths around R*. 



IV. ATOM-ION DYNAMICS 

Traps for individual ultracold atomic particles used in 
schemes for quantum information processing provide in 
most cases for the ability to manipulate the particles' 
motion via appropriate tuning of external trap parame- 
ters. This is in particular the case for optical lattices and 
Paul traps, where field polarizations and intensities can 
be changed to control the shape and position of the traps 
to a high degree of accuracy. Our proposal relies on these 
standard techniques, thereby introducing the innovative 
aspect of combining traps for ions and atoms. As al- 
ready discussed elsewhere [l^l, the physical mechanisms 
generating the traps for ions and atoms are different and 
lead, under appropriate conditions, to independent mi- 
croscopic traps, which can be modeled as follows. 

In this chapter dynamics will be described by intro- 
ducing a time dependent trap displacement d{t). Below 
a certain distance d ^ R* , the eigenenergies of the system 
start to depend on the spin state as well as on d itself, 
and positions of trap-induced resonances are determined 
by the internal state of both particles. In this way trap 
displacement can be used for spin-dependent control of 
the atom-ion system. 



A. Landau-Zener Theory 

The Landau-Zener formula gives a basic understanding 
of the atom-ion dynamics in the vicinity of trap-induced 
shape resonances. It describes a general two-level system 
whose eigenstates \ 'i>i) and |4'2) are coupled by some kind 
of interaction, and in which the two eigenlevels Ei and 
E2 form an avoided crossing when varying an external 
parameter. In our case this external parameter is the 
trap displacement d. The probability of a nonadiabatic 
passage of the crossing [l2| , 



exp 



(^i|g|^2)P 
h\ddEi2/dd\ 



(29) 



depends on the coupling matrix element, the velocity d of 
passage of the resonance and the relative slope dEi2/dd 



of the levels, where £'12 = £"1 — i?2- A fast passage of 
the avoided crossing (Pan ~ 1) results in a nonadiabatic 
evolution, preserving the shape of the wave function. 
At small velocities the resonance is passed adiabatically 
{Pna ~ 0), i.e. the system follows its eigenenergy curves. 
In the trapped atom-ion system, for certain trap sepa- 
rations, the energy of some molecular bound states be- 
comes equal to harmonic-oscillator energies, resulting in 
the avoided crossings. If such avoided crossing is passed 
adiabatically, then the initial harmonic-oscillator state 
with atom and ion located in their separated traps can 
evolve into a molecular bound state, where the atom and 
ion are trapped in a combination of the two external po- 
tentials. This process is reversible and we use it in Sec.lVl 
to realize an entangling two-qubit operation. 

In order to precisely predict the outcome of a col- 
lision process, in our simulations we have calculated 
time evolution numerically, using the Landau-Zener for- 
mula only as a guide to estimate the relevance of the 
avoided crossings for the transfer process. Since the en- 
ergy of the molecular bound states changes according to 
Emoi{d) ~ £:nioi(0)-h 1/2/2^2^2 (gee Fig. [3]), deeply bound 
states can cross vibrational states only at large d. In 
this case avoided crossings arc very weak, due to the fact 
that AE decays exponentially with the trap distance [l^l ■ 
Hence, the deeply bound states have no relevance for the 
dynamics and in our simulations we have included only 
shallow bound states that are closest to the dissociation 
threshold. 



B. Full dynamics in the single-channel model 

In the case of similar as and at, when the effective 
single-channel description is applicable, we describe the 
dynamics of the controlled atom-ion collision with the 
following time-dependent Hamiltonian: 



(30) 



We decompose the corresponding time-dependent wave- 
function in the basis of the eigenstates at d = 



|M/(t))=^c„,(t)|vl/ 



(0)\ 
nl h 



(31) 



in analogy to Eq. ([^5]) . Substituting this into Schrodinger 
equation, we obtain a set of coupled differential equation 
for the expansion coefficients c„;: 



nl 



^iuj'^d(t)D: 



n I 



(32) 



Here, D^'/ = {^^'i' l^l'^ni) dipole matrix clement, 

which in the context of the multichannel formalism is 
defined in Eq. ^7\ . We determine the radial part of the 



,(0) 



,(0) 
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single channel wavefunctions and the eigenenergies i?^^^ 
with the Numerov method [l^], using the effective short 
range phase as a boundary condition at minimal distance. 
From the wavefunctions, one can calculate the matrix 
elements . Inserting these into Eq. ([32]) . we are able 
to solve the equation for the coefficients numerically with 
standard routines. Thereby we verify in each case that 
our basis, limited by the maximal values of the I and n 
quantum numbers, is large enough and that the results 
do not change when increasing the basis. 

By comparing exact numerical dynamics with the re- 
sults predicted by the Landau-Zener theory, we have 
found, for example, an error of about 0.5% for the prob- 
ability of a fast, diabatic passage of the avoided crossing 
at d = 1 in the spectrum of Fig. [Sal with a speed of 1 
mm/s. Similar good agreement was observed in the adi- 
abatic limit, and only for intermediate speeds we found 
discrepancies of the order of 10%. This might be due to 
the complexity of the spectrum, which makes it impossi- 
ble to isolate an avoided crossing between two eigenstates 
from the influence of the rest of the eigenstates. Thus, the 
Landau-Zener theory is not applicable for quantum gate 
calculations performed for example in Sec. IV E II Also, 
faster processes lead to excitations to higher vibrational 
states that cannot be described by the Landau-Zener the- 
ory. 



V. QUANTUM GATE 

A. Qubit states 

In this Section we apply our model of the spin-state- 
dependent atom-ion collisions to construct a two-qubit 
controUed-phase gate. We encode qubit states in hyper- 
finc states of atom and ion. According to our previous 
notation, a given two-particle spin state is referred to as 
a channel a ~ {fifamf.mfjmi}. The total spin projec- 
tion m f = m f. + m is a conserved quantity during the 
collision. Therefore two states of unlike m/ cannot be 
coupled. In our case it is convenient to pick the compu- 
tational basis states 



|l>i,a = |/i.a=2,771/^^ =2),,Q 

according to Fig. [51 leading to the two-qubit states 

|00) = |/, = l,m/,=l,/a = l,m/„=l) = |1, 1,1,1) 
|01) = |/, = l,m/,=l,/„ = 2,r7v,=2) = |1,1,2,2) 

1 10) = |/, = 2,m/,=2,/„ = l,m/„=l) = |2,2,1,1) 

111) = |/, = 2,m/,=2,/„ = 2,m/„=2) = |2,2,2,2). 



(33) 



(34) 



Each of the two-qubit states is represented by a scatter- 
ing channel. The state |00) has m/ = 2 and is coupled to 
seven other channels, which have higher dissociation en- 
ergies and therefore they remain closed for |00) collisions. 
Thus the state |00) is stable with respect to spin-changing 



collisions. The channels |01) and |10), belonging to the 
ruf = S subspace, are coupled to each other and to two 
other channels, that are closed for both |01) and |10) 
collisions. There is no coupling for the state |11), since 
it is the only state in the = 4 subspace. In this 
way our choice of the qubit states minimizes the possi- 
bility of spin-changing collisions, and the only process 
that can lead to potential losses is the inelastic collision 
|10)-|01). 



i = I Atom (**^Rb) 



m-/a -2 -1 
/ = 2{ 

— 



|1> 



1^) 



i = I Ion C^'Ba+) 
m/, -2-10 1 2 
/ = 2{__ — — T) 



FIG. 5: (Color online) Specific choice of the qubit states 
out of the manifold of hypcrfine spin states of a *^Rb 
atom and a ^^^Ba + ion. 



B. Dynamics of four isolated channels 

Our multichannel theory describes collisions between 
atom and ion for general spin states of particles, in par- 
ticular for the four qubit states introduced in Eq. (|34p . 
For simplicity of the numerical calculations, here we do 
not perform the full multichannel dynamics, but rather 
we focus on the regime of applicability of the effective 
single-channel model, described in Sec. IIIIBI For ev- 
ery choice of parameters assumed in our calculations, we 
verify that the coupling to other spin states can be ne- 
glected. The total Hamiltonian including external de- 
grees of freedom, for the subspace corresponding to our 
computational basis, reads 

iJ = iJoo ® |00) (00| + i/oi ® |01) (01| 

-Hi/io® |10)(10| (8 |11)(11|. (35) 

We denote a qubit channel by \A) with A E 
{00,01,10,11}. Linear combinations of the computa- 
tional basis states form a general two-particle state j^*) = 
ayi|\I''^)|A), where denotes the quantum state 

of the atom-ion relative motion for the channel \ A). Ob- 
viously, the time evolution of 



-iEAt/hi-.j,A 



*-(i))|A) 



(36) 



is determined by the dynamics of the spatial part of the 
wave function |^'"^(i)), which we evaluate from Eq. (|32p . 
On the other hand, the phases due to the differences in 
threshold energies Ea in each of the channels, do not lead 
to state dependent dynamics, and can be eliminated by 
single qubit operations (see discussion in the next sub- 
section). 
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C. Phase gate process 

The two-qubit phase gate is represented by the trans- 
formation 



100) 


Interaction gi^oojQQ^ t/s^ 


|00), 


101) 


— > e^-^^lOl) — > 


101), 


|10) 


— > e'-^i^llO) — > 


|io), 


111) 


— > e'-^^lll) — > 


e^^ll), 



(37) 



performed on the computational basis states. The first 
step is the controUed interaction of atom an ion that leads 
to a specific phase for each two-qubit state. By applying 
the single-qubit transformation Us we can undo three of 
these phases and assign the total gate phase 

</> = </'oo + 4>ii - 001 - (j^w (38) 

to the |11) state [23|. If this phase equals tt, the phase 
gate, combined with single qubit gates, is a universal 
gate for quantum computation, as it is equivalent to a 
CNOT gate. It is possible to realize this phase gate 
scheme within our single-channel model, since the trans- 
formation of each two-qubit basis state can be treated 
separately. 

For our gate scheme, atom and ion are initially pre- 
pared in the motional ground state of their respective 
traps. The channel phases are gained by the control of 
the relative motion of atom and ion during the collision. 
Ideally we aim at obtaining back the motional ground 
state at the end of the gate process, so that the phase 
accumulated by relative motion is assigned to the qubit 
state. 

D. Gate fidelity 

Our definition of the fidelity is based on the overlap of 
the initial state of relative motion l^'j^J with the final 
state |\E''^(T)). In an ideal process these states are equal 
up to a state dependent phase. The fidelity needs to 
account for this phase. For one channel A and at zero 
temperature, according to [2^ we can define the fidelity 
Fa as follows 

Fa = 1 [1- |(*^J*^(T))|cos(7r- A0^)] , (39) 

where A(j)A ~ 4>a{T) — (t>'A is the difference between the 
desired channel phase 0^ and the phase (j)A{T) obtained 
by actual time evolution. In the following, we will assume 
that according to Eq. ([37]) the phases for the channels 
|00), |01) and |10) are undone perfectly due to the sin- 
gle qubit rotations leading to A0oo — A0io = A(/)oi = 0, 
while A(/)ii is nonzero. Hence, for channels |00), |01) and 
1 10) the fidelity is restricted only by the overlap between 
initial and final states, while for the state |11) we addi- 
tionally require that the total gate phase, computed from 
the single-channel phases with Equation psp . is = tt. 



We can further define the overall gate fidelity as 

f^gato = mm Fa, (40) 

since in our model the channels are decoupled (we neglect 
spin changing collisions). 

E. Adiabatic regime 

The adiabatic dynamics can be understood with the 
help of the correlation diagrams introduced in Sec. Ill CI 
Our gate scheme aims at an adiabatic transfer from an 
initial oscillator state to a molecular state ^'mob and 
back to the initial state. This is achieved by a variation 
of the trap distance across an appropriate avoided cross- 
ing, which we choose after investigating the correlation 
diagram. For example, the resonance at the trap distance 
d ~ 0.7 in Fig. I6al appears strong enough and we use it 
in numerical calculations in the following. During the 
transfer process each logical basis state acquires a differ- 
ent phase, since the energies of molecular states depend 
on the channel (see Fig. [Gb]) . The phase accumulated for 
each channel in an adiabatic transfer process is given by 
the integral 

0pot = -^/'°'^i?^Wdt, (41) 

where is the energy of the adiabatic eigenstate de- 
picted as a function of d in Fig. [5] (thick curve). In the 
adiabatic regime we avoid excitations to higher vibra- 
tional states by keeping the velocity d small, compared 
with the characteristic velocity of the harmonic motion 
in the trap: d ^ yjhui/ ji. On the other hand we choose 
the velocity across weaker resonances at larger distances 
high enough to pass them diabatically, as seen in Fig. [B] 
We want to find a particular function d{t), which, if ap- 
plied on the trapped atom-ion system, results in a de- 
sired total gate phase, while ensuring diabatic passage of 
the weak resonances as well as adiabatic passage of the 
strong resonance. Since the total phase depends on the 
difference between single-channel phases (see Eq. 
the gate speed, in fact, is determined by the differences 
in the potential energy curves of unlike channels. 

For the sake of concreteness we assume specific values 
of the singlet and triplet scattering lengths, in such a way 
that our single-channel effective model is applicable. For 
our calculations we choose as = 0.90i?* and at = 0.95i?*. 
According to the procedure described in Sec. lIII Bl the es- 
timate of the error introduced by the model is 2 x 10^"^. 
For singlet and triplet scattering lengths that differ by 
more than 10% the channel mixing becomes already sig- 
nificant and does not allow for a single-channel descrip- 
tion. 

Actually, the singlet and triplet scattering lengths are 
uniquely determined by the specific choice of the atom- 
ion system we describe. So far these parameters have not 
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been measured experimentally, for any atom-ion system. 
However, as soon as the accurate values of Os and at 
arc determined, one can repeat the calculations with the 
physically correct parameters, which may require going 
beyond the single-channel model, and including the full 
multichannel dynamics according to Sec. IIIII 

Assuming the single-channel effective model we first 
compute a correlation diagram for each of the channels. 
This is done by diagonalizing the Hamiltonian in the ba- 
sis of eigenstates evaluated at ci = 0. The result is de- 
picted in Fig. [Sa] for the |11) channel. The diagrams 
show small differences in the molecular states at small 
distances (Fig. [Gb]). since the energy of molecular states 
depends on the atom- ion spin configuration (the qubits). 




0.4 0.6 0.8 1 1.2 1.4 1.6 0.59 0.6 0.61 0.62 0.63 



d/B' d/R' 

(a) (b) 

FIG. 6: (Color online) Correlation diagrams for = 0.90 
and at ~ 0.95, where for each of the qubit pairs we 
subtract the threshold energy of the corresponding chan- 
nel. The energy curve followed in the adiabatic process 
is marked with a thick, red (dark gray) line in the left 
Figure (a). We only show the complete diagram for the 
|ll)-channel. Small energy differences of the channels 
around d = dmin can be seen in the close-up (b). These 
differences are the basis of our proposal for realizing an 
atom-ion phase gate. 



the other hand, much larger velocities would cause un- 
wanted motional excitations in the trap. At d = 0.95i?*, 
the velocity is decreased to = O.IR* /{h/ E*) in order 
to adiabatically convert the trap state into a molecular 
bound state using a stronger resonance. The curve is 
followed down to some minimal distance dmin- Then, 
the reversed pulse brings the system to the initial trap 
separation. Fig. [7a| shows the complete d{t) function. 
It is known that sharp kinks can cause motional ex- 
citations, therefore we use a smooth function d{t) = 
d -f 1/2 (d'l + 4)^ ± ^/(d'~~d!^fFTC^ to change be- 
tween two slopes d'i 2 ~ 0.5 and Q.IR* / {h/ E*). Here, 

the +(— ) sign is used for incrcasing(decrcasing) slope, d 
is an offset and C is a parameter adjusting the curvature 
at the kink. Around the turning point at dmin d{t) ~ . 

We can now compute the solution of the Schrodinger 
equation numerically at a given time by solving Eq. (|32p 
with standard numerical routines . This yields the gate 
phase which, for example, can be adjusted by a variation 
of dmin- This phase is in fact a phase difference accumu- 
lated due to the energy splitting 

= S,™ 1 + Ell^^ - E^^^ - E'^^^ (42) 

shown in Fig. [8] function of trap separation; the 
larger Ai5, the faster a phase difference is reached. Thus, 
decreasing dmin increases the phase, as seen in Fig. [Tb] 
We find a gate phase ot cf) = I.OOQtt at dmin = 0.591i?*. 
The corresponding gate fidelity is i^gatc = 0.994 according 
to Eq. (gOl)- The process takes the time T = 9.14;i/£;*, 
which equals T = 1.31ms for our choice of the trapping 
frequency w = 27r x 30kHz and the masses and hyperfine 
structure of the *^Rb atom and the ^•^^Ba+ ion. We show 
the population of the instantaneous eigenstates in Fig. O 
At c?max the trap ground state is labeled with n = \. The 
system is initialized in this state. At half gate time the 
quantum state changed to n = 4, which is the molecular 
state at dmin- Finally, at the end of the gate process, we 
almost perfectly obtain back the initial state. 



1. Numerical simulation of an adiabatic phase gate 

In our simulation of the gate process the initial and 
final trap separations coincide: d(0) = d{T) = (i,„ax- 
We assume that initially the atom and ion are each in 
the ground-state of its own trap. The distance dmax is 
determined in such a way that there are no bound states 
in the vicinity of dmax that would influence the harmonic- 
oscillator ground-state. 

The controlled time evolution of our atom-ion system 
requires the appropriate adjustment of the distance d{t) 
between the two trapping potentials as a function of time. 
The slope of this function is essential for the result. 

In order to follow the energy curve depicted in Fig.[6a[ 
we construct a specific function d{t). We start at dmax 
with an initial velocity d'l = 0.5R* / {H/ E*), which is large 
enough to traverse weaker resonances diabatically. On 



F. Fast gate using optimal control 

Quantum optimal control techniques are a powerful 
tool that allows to increase the fidelity of a time evolu- 
tion process by finding an appropriate pulse shape for 
some external control parameter. The outcome of the 
controlled collision of an atom and an ion is very sensi- 
tive to the particular shape of the time dependence in 
the trap distance d{t). It is hard to manually design a 
specific function d{t) that leads to a phase gate fidelity 
very close to unity. 

The most significant problem in our specific example 
calculation results from a relatively weak resonance at 
do ~ IR* , which we want to pass in a non-adiabatic way. 
We could not find an optimal constant slope that brings 
us from the trap state at d > dp to the trap state with d < 
do across the resonance without losses. Large velocities 
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(a) (b) t(m,) 



FIG. 7: (Color online) (a): The function d{t) starts 
at dmax = l-4i?* with a slope of d[ = 0.5R*/{h/E*) 
large enough to diabatically pass weaker resonances. The 
velocity is changed to around d = 0.95R* to d'2 = 
Q.IR* / [h/ E*) in order to ensure an adiabatic traversal of 
a stronger resonance. The system is brought to the ini- 
tial distance with the reversed pulse and the kinks are 
smoothed to avoid motional excitations. The charac- 
teristic unit of speed is R*/{h/E*) = 2.05mm/s. (b): 
gate phase as a function of d^i^ using the described d{t) 
shape. We find that with d^-^ = 0.591 a gate phase of 
(p = I.OOOtt is reached. 
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d/R* 

FIG. 8: (Color online) Energy difference as a function of 
trap separation. For d < 0.3R* the function is practically 
constant, which, for our phase gate process, means that 
bringing the traps closer than dy^in ^ 0.3i?* does not lead 
to a significant advantage. 



FIG. 9: (Color online) Population of the most impor- 
tant adiabatic eigenstates for the qubit channel |11) dur- 
ing the adiabatic gate process (the remaining channels 
show similar behavior). States are labeled according to 
their energy in the correlation diagram in ascending or- 
der. The initial state at t = is the relative motion 
ground state in the trapping potential for d = (imax- The 
label of this state is set to n = 1. The molecular state 
marked in Fig. [6a| then bears the label n = 4. Around 
a crossing, the labels of molecular/trap states switch be- 
cause their energies change order. Our goal is to follow 
the thick red (dark gray) curve in Fig. [Bal We observe 
that indeed two resonances are passed diabatically. At 
t ~ T/2 the state n = 4 is reached with relatively high fi- 
delity, which means that a molecular ion is formed here. 
At d{T) = dinax we regain the initial state again (the 
state dependent phase of the final state is not depicted 
here). To better show the features of the curves, the time 
axis is squeezed between ti = 0.12 ms and t2 = 1.18 ms, 
where the velocity is lower (see Fig. [7|). 



with an initial guess for the control function d(t), which 
in general does not yield a satisfactory fidelity. The fi- 
delity of the process is increased in every iteration step 
by updating d{t). We divide the time axis in small time 
steps dt. At each time step we evolve the wave function 
forward in time using the Crank- Nicholson scheme [2^ . 
The update of d{t) is successively performed in every time 
step. 



lead to excitations of energetically higher states, while 
the consequence of smaller velocities is a non-negligible 
population of the molecular state that crosses the trap 
state. This population is not fully recovered on the way 
back. In Sec. IV E II we nevertheless found a process that 
yields a gate fidehty of Fgato = 0.994. 

With optimal control we can not only find a pulse 
shape that produces a satisfactory gate fidelity, but also 
go beyond the adiabatic regime and reduce the gate time. 
By applying larger velocities, we allow for excitations to 
higher energy levels. Making use of interference effects, 
an appropriate d{t) pulse shape can undo these excita- 
tions in the final state of the process. 

This optimal pulse shape is found here with an itera- 
tive optimization algorithm called intermediate feedback 
control, which is for example introduced in (isj . We start 



1. Enhancing the fidelity of the adiabatic gate 

The adiabatic gate process of Sec. IV E II has the fi- 
delity Fgato = 0.994. With only three iterations of the 
optimal control algorithm we can enhance this fidelity to 
Fgatc = 1 — 7 X 10'* , which even exceeds the validity of 
the underlying single-channel model. The gate phase is 
improved to 1.00267r. The optimized function d(t) shows 
small scale variations ('wiggles') that arc a typical feature 
of the used optimization algorithm. These wiggles have 
the amplitude ~ 0.002i?* « 0.6nm and happen on a time 
scale of the order of 10 ps. This amplitude is smaller 
than the uncertainty of the ion trap center position in 
up-to-date experimental realizations. 
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2. Fast Gate Scheme 

It is desirable to reduce the gate time to a minimum. 
In our case this minimum is given by the time that is at 
least required for accumulating the gate phase = tt. We 
profit from the energy differences of the molecular states 
for the different channels. The differences are largest at 
small distance dmin, but below d « 0.3i?* they are prac- 
tically constant. The fastest possible gate needs to effec- 
tively transport the atom-ion relative wavefunction into 
a molecular state, where the gate phase is accumulated 
during a certain time. The gate ends with a reversed 
pulse and brings the system to the initial trap ground 
state, while the logical phase is preserved. 

The optimal control algorithm can be used to build 
such a gate process. As a first step we design d{t) 
changing the quantum state of atom and ion from the 
trap ground state at dmax = lAR* to the molecu- 
lar state at dmin = 0-3R*, for all channels during the 
transport time itrans- The optimization objective J ~ 
X;^2Re{(^'^(itrans)|*nioi>} ^^r this step aims at max- 
imizing the overlap of the time-evolved state with the 
desired molecular state. The second step is a stationary 
evolution at dmin, where the main part of the differen- 
tial phase is accumulated. Subsequently we perform the 
reverse of the initial pulse. The combined d{t) pulse is 
shown in Fig. llOal 

The total gate time is T = 346ps. With respect to 
the adiabatic case this means a reduction by a factor 
of 4. Since our scheme makes use of the channel energy 
differences in the molecular state at d — dmin ~ 0.3i?*, 
we can estimate a quantum speed limit for this process 
Tiimit = ttH/AE, with AE from Eq. (gl]), which gives the 
minimal gate time, neglecting transport durations and 
infidelities of the molecular state's population. In our 
example this limit is Tumit ~ 250ps. Our gate process 
time lies very close to this value, considering an overall 
transport time of 2 * ttrans = 158ps. We note that a 
part of the phase is accumulated during the transport 
phase, since we already enter the regime where the energy 
differences become significant. 

Further reducing the transport time may be possible, 
but the optimization algorithm we used stopped con- 
verging in reasonable time for larger transport veloci- 
ties. However, we already could significantly accelerate 
the process utilizing optimal control techniques. 



3. Perspectives for further improvement 

Certainly the gate speed would improve if the differ- 
ence between singlet and triplet scattering lengths was 
larger than assumed in our example calculation. In gen- 
eral case one can use our multichannel formalism for an 
accurate description of the dynamics, beyond the single- 
channel approximation. Wc point out that the essential 
parameters and at are still unknown and they need to 
be measured in experiments in order to do realistic calcu- 



lations for specific systems. However, our work demon- 
strates the feasibility of an ion-atom quantum gate even 
based on the simplified scheme we assume for calcula- 
tions. 

Further possibilities occur for very different values as 
and at. In this case especially the |10) and |01) states 
are coupled strongly. In this case optimal control mech- 
anisms can be used to suppress effects of spin-changing 
collisions in the final state, in order to realize a phase 
gate. The coupling of |01) and |10) could also be effec- 
tively used for a SWAP gate, which inverts the popula- 
tions of these two channels — or for its square root, which 
in combination with single-qubit rotations constitutes an 
alternative universal set of gates for quantum computa- 
tion. 



VI. CONCLUSIONS AND OUTLOOK 

In this work we analyzed the spin-state-dependent in- 
teraction between a single atom and a single ion guided 
by external trapping potentials. We applied our insight 
on this system in order to realize a two-qubit quantum 
gate process and thereby provide the basic ingredients 
for quantum computation with atoms and ions combined 
in one setup. This work was motivated by recent exper- 
imental possibilities combining magneto-optical traps or 
optical lattices for atoms, and RF traps for ions. These 
experiments are currently established in several groups 
worldwide 

Wc started our description of controlled interaction 
of an atom and an ion by formulating a multichannel 
quantum-defect theory for trapped particles, analogous 
to the free-space case discussed in Q. This step simpli- 
fies the description of atom-ion collisions as it does not 
require a detailed knowledge of the molecular potentials 
at short range. Experiments measuring the positions of 
Feshbach resonances can determine the essential two pa- 
rameters for our theory — the singlet and triplet scatter- 
ing lengths. Since these experiments have not yet been 
performed, here, we have focused on the case of simi- 
lar singlet and triplet scattering lengths given in order 
of magnitude by R*, while the general case of different 
scattering lengths have been discussed only qualitatively. 

We were able to reduce the multichannel formalism to 
an effective single-channel model that singles out a spe- 
cific spin state of atom and ion. This model is found to 
be accurate for similar values of the singlet and triplet 
scattering lengths. In calculations we have assumed 
as = 0.9R* = 4989ao and at = 0.95i?* = 5266ao for 
the singlet and triplet scattering lengths, respectively. In 
contrast, opposite signs of scattering lengths exclude a 
single-channel description. We estimated the error intro- 
duced by our specific single-channel model to 2 x 10^^, 
which is due to a mixing of channels in the eigenstates 
of the system. Taking even closer values of singlet and 
triplet scattering lengths lead to a better applicability of 
the single-channel description. 
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FIG. 10: (Color online) (a): Optimized c?(t) funetion and 
initial guess (dashed) for the fast gate process. The ini- 
tial state is the trap ground state at dmax- The popula- 
tions of the adiabatic eigenstates are changed by the opti- 
mization process. Applying the initial guess pulse higher 
vibrational and also molecular states become populated 
(b) . These excitations are prevented or undone using the 
optimized pulse (c). In this case, at dmin the desired 
molecular state n = 4 is reached almost perfectly. For 
simplicity we only show the plots for the qubit channel 
1 11) here, the situation is very similar for the other chan- 
nels. Note the different time scales used for the transport 
and phase-accumulation sequences respectively. 



Where applicable, our effective single-channel model 
can be implemented in calculations in the context of ul- 
tracold chemistry as well as ultracold scattering physics. 
A single-channel description assigning quantum-defect 
parameters to each channel separately was already dis- 
cussed in [12}. However, the spin state was not included 
in previous research. In the present approach, starting 
from the fundamental parameters Og and at of the multi- 
channel formalism, we derive the quantum-defect param- 
eters of each isolated channel consistently. We take the 
channel coupling into account and we estimate the error 
introduced by assuming isolated channels. Therefore, we 



can apply the model to quantum computation schemes 
that store qubits in internal spin states of atom and ion. 

A remarkable feature of the system is trap-induced 
shape resonances that couple molecular bound states 
to unbound trap states. Quasistatic eigenenergy curves 
show these resonances as avoided crossings. They can be 
used to form ultracold trapped molecular complexes and 
thereby allow full control of cold chemical reactions. 

Trap-induced resonances form a basis for our idea of 
the phase gate process as well. Initially an atom and 
an ion are prepared in the trap vibrational ground state. 
We realize a qubit-dependent two-particle phase via con- 
trolling the external degrees of freedom. By bringing 
the traps close together we let the particles interact and 
finally separate them, obtaining the motional ground 
state again. In doing so we cross weaker resonances di- 
abatically (remaining in a trap state) and then follow a 
stronger resonance adiabatically into a molecular bound 
state, where a two-qubit phase is accumulated. Since 
the positions of the resonances are different for each spin 
combination, the accumulated phase is different for each 
qubit channel and we are able to control the trap dis- 
tance in such a way that a two-qubit phase gate is real- 
ized. This phase gate, in combination with single qubit 
rotations, is a universal gate for quantum computation. 

We performed numerical simulations of the controlled 
collision specifically for a ^^^Ba"*" ion interacting with a 
^''Rb atom, each guided by a spherically symmetric har- 
monic trap with uji^a = 27r x SOkHz. We have chosen 
specific hyperfine-qubit states to obtain the four qubit- 
channels 00,01,10,11. In this framework we developed 
a two-qubit phase gate process entangling atom and ion. 
We thereby showed that trap-induced resonances can be 
used to control the interaction of atom and ion. The er- 
ror for our gate process is 1 x 10^'^ and in this case the 
gate time is 1.3 ms. Using optimal control techniques we 
were able to accelerate the process to 346]is. In future 
work we plan to decrease the gate time by using higher 
trapping frequencies allowing faster transport. 

Our choice of very similar scattering lengths allowed 
to perform single-channel calculations, but on the other 
hand it limits the gate time, as the energy splitting be- 
tween the channels is rather small. A more general de- 
scription can be done in the framework of our multi- 
channel formalism, allowing for arbitrary combinations 
of singlet and triplet scattering lengths. It is possible 
that the actual values of the scattering lengths are in fact 
very different, which would require a more complicated 
multichannel computation, but also possibly allows much 
faster quantum gates. However, already in the regime 
we considered, a gate time below a millisecond is demon- 
strated. 

In our model, we assumed harmonic trapping poten- 
tials, and the oscillator frequencies were equal for atom 
and ion. Among the advantages of ions for quantum 
computation is the existence of much tighter trapping 
potentials than for atoms. The basic ideas developed in 
this paper are expected to be applicable to such more 
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general situations with different trapping frequencies. 
A generalization is highly desirable, but will lead to a 
more complicated theoretical treatment, as for example 
center-of-mass motion becomes coupled to the relative 
motion. Very elongated cigar-shaped traps have already 
been treated in previous works. One of our goals is the 
consideration of particular experimental realizations in 
order to describe them with our theory, compare the re- 
sults or suggest directions of experimental research. 

In the present work we did not use external mag- 
netic fields in order to manipulate the interaction. Mag- 
netic induced Feshbach resonances have been applied 
very successfully to engineer neutral atom collisions. Fu- 
ture investigations will include magnetic fields to control 
the atom-ion interaction even more efficiently and possi- 
bly combine trap-induced resonances and Feshbach res- 
onances for this purpose. Our work can be seen as a 
principle investigation of a new, interesting physical sys- 



tem and can be extended in many directions. 
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